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Abstract 

The interaction of electrically charged particles in a dilute gas of pointlike 

magnetic dipoles is studied. We show that the interaction potential at small 

distances has a linear piece due to overlap of the dipole clouds gathered near 

electric sources. At large distances the potential becomes of the Coulomb type 

with nonperturbatively renormalized charge of the test particle. The physical 

applications of these results are discussed. 
14.80.Hv, 11.15.Tk 
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I. INTRODUCTION 



We study properties of electric charges in the plasma of the Abelian magnetic monopoles 
in three dimensional Euclidean spacetime. Physically, the magnetic dipoles may arise as 
the monopole-antimonopole bound states in gauge theories like, for example, the Georgi- 
Glashow model. This model possesses the topologically stable classical solution called 
't Hooft-Polyakov monopole [[I]] which carry a unit of magnetic charge (the instanton num- 
ber in three dimensions). The 't Hooft-Polyakov monopole consists of a compact core and 
longrange gauge fields associated with unbroken Abelian subgroup. 

In the weak coupling regime the vacuum of the Georgi-Glashow model is filled up with a 
dilute monopole-antimonopole plasma. Due to a longrange nature of the Abelian fields the 
behavior of the plasma is essentially the same as that of the Coulomb plasma of the Abelian 
charges. Polyakov has shown that in this plasma test particles with opposite electric 
charges experience confining forces at large separations due to formation of a stringlike 
object between the charge and the anti-charge. The string has a finite thickness of the order 
of the plasma correlation length (i.e., the inverse Debye mass) and a finite energy per unit of 
the string length ("string tension"). Thus the potential between the test particles is linear 
at large distances 0. 

Small distance behavior of the inter-particle potential in the monopole gas is also inter- 
esting. According to Ref. || at distances R much smaller than the correlation length of 
the plasma the potential contains a nonperturbative piece proportional to R a , a ~ 0.6 in 
addition to a perturbative contribution due to the one photon exchange. The former is a 
consequence of the string formation at large distances. A nontrivial short distance potential 
may have many physically interesting consequences, see, e.g. Ref. ||,]5] for a discussion in 
context of QCD and other theories. 

At a sufficiently high temperature the Coulomb plasma was shown to undergo the Bere- 
zinsky-Kosterlitz-Thouless phase transition ||[7| . In the high temperature phase the charges 
form neutral bound states 0HH which obey nonzero dipole moments. Thus at the phase 
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transition the Coulomb plasma of charges transforms into a gas of the dipoles in which 



the Debye screening is absent pO|,p[]. The field of the magnetic dipole is too weak at large 
distances to be able to induce a nonzero string tension between the electric charges. As a 
consequence of this fact, at large separations between the electric charges the string structure 
is destroyed in the high temperature phase and the potential becomes of the Coulomb type. 
A detailed study of the temperature phase transition in the three dimensional Georgi- 
Glashow model can be found in Ref. ||. The physics of dipoles in gauge theories has also 



been discussed in Ref. 11 



Despite of the absence of the longrange confinement the physics of the electric charges 
in the dipole gas is still interesting. In this paper we show that the short range potential 
between the charges has a linearly rising piece due to interaction of the electric charges 
with the magnetic dipoles. For simplicity we consider a dilute gas of the pointlike dipoles 
characterized by a (fluctuating in general case) dipole moment fl and the spacetime position 
x. The temperature effects are not studied in the present publication. The structure of the 
paper is as follows. The path integral formulation of the dilute dipole gas is presented in 
Sec. [TT| The interaction of the magnetic monopole gas with the electric charges is discussed 
in Sec. fTJ and interaction energy of the static charge-anticharge pair is studied in Sec. [TV] . 
The generalization of the present approach to the fluctuating dipole moments is given in 
Sec. [VJ Our conclusions are summarized in the last section. 



II. PATH INTEGRAL FORMULATION 

The magnetic dipole is a magnetically neutral localized pair of monopole and anti- 
monopole separated by the distance r. The magnetic moment of the dipole is ft — g m • r 
where g m is the magnetic charge of the constituent monopole and r is the distance between 
the constituents. If the typical distance between the dipoles is much larger than the dipole 
size r then the dipoles may be treated as pointlike particles. This condition can be written 
as follows: 
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p3r«l, (2.1) 

where p is the mean density of the dipole gas. In this case the dipole is characterized solely 
by the position x and the dipole moment p while the "dipole size" r becomes an internal 
characteristic of the dipole. The action of two interacting pointlike dipoles with magnetic 
moments p a and p b located at positions x a and x b , respectively, is given by the formula: 

V(p a , p b ; x a , x b ) = (p a ■ d) (p b ■ d)D {3D) (x a , x b ) , 

where D^ 3 ^(x) = (47r|x|) _1 is the propagator for a scalar massless particle in three dimen- 
sions. 

The statistical sum of the dilute dipole gas can be written as follows: 

°° C N f f f f [ 1 N - ] 

2 = J d3x i J '"J d 3 XN J d 3 p N expl -- V(p a ,p b ;x a ,x b ) \ , (2.2) 

N ° a^fe 

where ( is the fugacity parameter Q. As in the case of the monopole gas the dipole 
fugacity is a nonperturbative quantity since ( ~ e~ s °, where S ~ g~ 2 is the action of a 
single dipole and g e = 2Tr/g m is the fundamental electric charge in the theory. 

The integration over the dipole moment / d 3 p in eq. (|2.2|) is given by the integration 
over direction n of the dipole moment p a = p a ■ n a (no sum is implemented) and over its 
absolute value, p a , weighted with a distribution function. For the sake of simplicity we 
fix the absolute value of all magnetic moments p a = p. The case of the fluctuating p is 
considered in Sec. |V|. 

Introducing the scalar field \ we rewrite expression ( |2.2| ) as follows: 

exp{- K / D Xexp{-i J d 3 x0x) 2 

— CO 

N 

+ipJ2(™a ■ d)x(x a )} , (2.3) 



a ' b y^ -oo 

N 



0=1 



1 Note that in 2 + 1 theory the dimensionalities are: [g^] = mass 1 , [p] = mass 2 and [£] = mass 3 . 
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where (a,b) is the scalar product of the vectors a and b. Substituting eq.(|2.3|) into eq.(pi 
and summing over all dipole positions, we get the following representation of the dipole 
partition function: 



Z oc J Dx ex p{~2 J d3x $X) } ' Z I J d3x J^ n ex P{^( ? 

— n 2 =l 

= / Dx exp{-y d 3 x \x0xf - C J d 3 n exp{i^(n ■ d)x{x)}]} . (2.4) 



N 

(n • d)x{x)} 

,7-: 'I 



After integration over the vector n we obtain (up to an inessential multiplicative factor): 

Z = j D X exp{- J d 3 xC[\d X \]}, (2.5) 



where \dx\ = v <9x > an d 

A/) = ^ 2 -4vrC^W), «(/) = 5^. (2.6) 

We treat the model ( j2.5| , [2.6| ) as an effective theory. Shifting the field / = \dx\, f —> f / 

- _i 
\i = f and rescaling the space coordinates, x — > (47r^) 3 x we immediately determine that 

the dynamics of the model is controlled by the dimensionless constant 

A = 4vrC/i 2 , (2.7) 
and the typical field fluctuation is of the order of f 2 ~ A. Thus at small values of the 



parameter A lagrangian Q2.6|) can be expanded into series of f 2 and the interaction terms 
f 2n with n > 2 can be treated as small perturbations. 

The vacuum expectation value, p =< p(x) >, of the dipole density, 

p(x) = 5]5 (3) (a;-a; i ), 



can be rewritten in representation 

p=iw({ ^mm )=4<l+0{x)): (28) 

\dx{x)\ 

where the last equality is written for small couplings A. In this regime the coupling A is 
proportional to the density of the dipoles, A = p p 2 + 0(A 2 ), and therefore the condition 
A 1 can be interpreted as a requirement for the density of the dipole moments to be small. 



III. ELECTRIC CHARGES IN MAGNETIC DIPOLE GAS 



An infinitely heavy test particle carrying the electric charge qg e in introduced in the 
vacuum of the theory with the help of the Wilson loop operator, W q (C). Here contour C is 
the trajectory of the particle. Let us first consider the contribution of the Abelian monopoles 
to the Wilson loop: 

W q {C) = exp{z J d 3 xp mon (x)r] C (x)} , (3.1) 

where p m0 n(x) = 9m J2a m a^(x — x a ), m a G % is magnetic charge density. The function r)c 
is defined as follows: 

rf( x ) = n J d 2 a i:j (y) e ijk d k D(x - y) , 

where the integration is taken over an arbitrary surface Sc spanned on the contour C. Note 
that the value of the Wilson loop does not depend on the shape of the surface E<?. 

Now suppose that all monopoles and anti-monopoles appear always as dipole pairs of 
the small size r. In the pointlike dipole approximation eq. ( |3.1| ) can be rewritten as follows^: 

W q (C) =exp{trq J d 3 x (^(x), d)7] C (x)} , 

= exp{irq J d 3 x (^(x), 7 C (a;))} , (3.2) 

where 7 C = dr] c , and ffa\x) = J2i P-i^ 3 \x — x«) is the density of the dipole moments. Since 
the charge of the Wilson loop q always enters eq.( |3.2j ) in the combination qr we consider 
below the Wilson loop for the unit charge, W(C) = W\{C), and then make the rescaling 
r — > qr in the final result to describe the potential for q ^ 1. 

To study the static potential between charges we consider the flat Wilson loop operator 
defined in the plane (xi,X2), see Fig. [I]. The operator inserts the static charge and anti- 
charge pair separated by the distance R. The particle trajectory consists of the two straight 

2 The function rjc lias the jump r)c —> f]c + 2tt at the position of the surface Sc. This jump is 
inessential for our considerations due to 27r-periodicity of the exponential function of an imaginary 
argument. 
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lines, C = L\ + L2 and the surface £e is located in the plane (x\,X2)- The corresponding 
function rj c is: 

c, \ ( x \ + R / 2 \ fX\ - R/2\ 
ri (xj = arctan — arctan , 

v £3 ' v x 3 ' 

and the function 7 C is defined as follows: 

/0 -* 
d 2 a ij (y)e ijk —d x D (x-y) (3.3) 

TsingO*') Treg("^) i 

iLzM = * ^3 6{x 3 ) 0(2x! + i2) 0(2x! - i2) ; (3.4) 



f ri ' (x) =AR 8 xi z 3 , , R 2 + 4(zij - ^ 



x 



-l 

2 



~(i? 2 + 4(x 2 + :r 2 )) -16i? 2 x 2 " 

where G(x) = 1 if x > and Q(x) = otherwise. 

Substituting eq. ( |3.2|) into eq.( |2.2| ) and performing the transformations presented in Sec. 
Il| we represent the quantum average for the Wilson loop in the dilute gas as follows: 

< W(C) > gas = ±Jd X exp{-| d 3 x ^-0 X f 
+AnC(l-v(\fid X + rf\))}}, 

where the expression in the square brackets is normalized to zero at x = and the function 
v is defined in eq.( ^6|) . 



Physically, the magnetic dipole gas affects the potential between oppositely charged elec- 
tric particles because of overlapping of the dipole clouds gathered near the electric sources. 
Indeed, from the point of view of the three dimensional physics the static test charge tra- 
jectories C = Li + L 2 , Fig. |I|, may be considered as an electric current running along the 
contour C. This current induces a magnetic field which circumvents each of the paths L\ 
and L 2 . Since the magnetic dipole lowers its energy in the magnetic field, the density of the 
magnetic dipoles should increase towards the position of the test electric charges. The in- 
teraction energy of the dipole clouds gathered near the oppositely charged particles depends 
on the inter-particle separation. Below we show that at small charge separations this energy 
is a linear function of the distance between the charges. 



IV. POTENTIAL OF STATIC CHARGE ANTICHARGE PAIR 



The static interaction potential of particles with electric charges ±gg e is a sum of the 
perturbative contribution from the one photon exchange and the contribution from the 
dipole gas, respectively: 

V q {R) = -^f D {2D) (R) + E*»(R) , 
E^(R) = - lim i < W q (C RxT ) > , (4.1) 

1 — >oo ± 

where D^d){R) — — (2tt) _1 log(mR) is the two dimensional propagator for a scalar massless 
particle, m is the massive parameter and Cr x t stands for the rectangular R x T trajectory 
of the test particle. In the limit T — > oo the contour Cr x t transforms into the contour 
C = L x + L 2 , Fig. |. 

We evaluate the dipole gas contribution E& as classically minimizing the functional (we 
put q=l), 

l tZ \2 , a sfi sin(\ /I dx + r ^ \)~ 



S{ X ) = Jd 3 x [-(9 X ) 2 + 4vrC(l 



with respect to the field x i n the limit A € 1. To this end we note that the singular part 
( 3.4|) of the source (|3.3|) does not contribute to the functional S(x)- This statement can be 
checked directly either by regularizing the singular piece of the source 7 at a finite R or by 
considering the limit R — >• 00 in which only the singular part 7f ing is nonzero at finite X\ 
and £3. Applying rescaling arguments similar to those mentioned after eq.( [2.6|) we conclude 
in the limit A 1 the classical solution is suppressed as \dxcA = 0(X) and therefore in the 
leading order in A the classical energy in the leading order is given by the formula: 

Ef s (R) =4ttC Jdx! Jdx 3 (l-v(rT c (x,R)))] , (4.2) 
r c (s, R) = 4R [(r 2 + A(x 2 3 + x\)f - \QR 2 x]\ " 3 . (4.3) 

Unfortunately the explicit analytical integration is impossible due to the complicated form 
of the function T, eq.fl4.3p. However, if the distance R between the test particles is much 
smaller than the "dipole size" r then the energy E C \(R) can be expanded in series of R/r: 
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E^ S (R) = 7T 3 (rR(l + 0(j)), R < r . (4.4) 

Thus, the energy of the charge-anticharge pair is the linearly rising function of the inter- 
particle distance 

Note that in realistic theories the magnetic dipoles are realized as monopole-antimono- 
pole bound states with a nonzero physical size r. If the distance between test particles R is 
much smaller than the size of the such dipole, R <C r, then the test particles become affected 
by the monopole constituents and the pointlike approximation to the dipole is no more valid. 
However as we will see in the next Section the linear potential extends over R ~ (a few) r 
and therefore we believe that our results derived in the pointlike dipole approximation may 
also be applicable to real systems. 

According to eq.( [4.4|) the coefficient a q in front of the linear term in the interaction 
potential of the particles carrying the electric charges ±qg e is: 

a q = qa u en = n 3 (r = — pr , (4.5) 

where p is the dipole density, eq.( |2.8| ). This formula is given in the leading order of A. To 
derive eq.( [4.5| ) we have used rescaling r —>■ qr according to the discussion after eq. (|3.2| ). 

The coefficient a q is proportional to the electric charge of the test particle. Thus, each 
elementary electric "flux" coming from the charge is not interacting with the other fluxes 
and, as a result, the total coefficient, a q , is a sum of the elementary coefficients <i\ of q 
individual fluxes. This is due to the fact that leading contribution to the formation of the 
dipole clouds is given by the external magnetic field induced by the test particles. The dipole 
dynamics would affect the dipole clouds in next to the leading order terms in expansion over 
the A parameter. Since the density of the dipole clouds is a function of the test particle 
charge, then the dipole density-density interactions in higher order in A would lead to q n , 
n > 1 corrections to the linearity coefficient a q . Thus it is natural to suggest that in general 
case the proportionality of the linear term to the electric charge of the test particles is lost 
even on the classical level. Note that the proportionality of the linear term coefficient a 
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to the number of elementary fluxes appears in the so called Bogomol'ny limit [12 of the 



Abelian Higgs model for the Abrikosov-Nielsen-Olesen strings |I3[. However, in the latter 
case the proportionality of the string tension to the number of fluxes is exact on the classical 
level. 

At large distances, R 3> r, the dipole gas contribution to the interaction energy ( |4.2| , [4.3D 
between the oppositely charged electric particles grows logarithmically: 

ET{R) = ^Cr 2 log(-) [l + oQ] , J2>r, (4.6) 



in accordance with absence of the charge screening in the dipole gases [|K]J§. Shifting the 
dipole size r — > qr in eq.( |4.6|) we conclude that the dependence of the nonperturbative 
potential ( |4.6|) on charge of the test particles, q, and on the distance between them, R, is 
essentially the same as for the perturbative one photon exchange, eq.( [4.1| ). Thus the dipole 
gas nonperturbatively renormalizes the coupling constant g e at large distances and the full 
potential ( |4.1| ) has the following behavior: 

2 2 

VJR) = log(R) + const. , 

4:71 

e=l + \x + 0(X 2 ), i?>r, 

where e is the dielectric constant. Note that both the appearance of the linear potential at 
small distances and the renormalization of the electric charge at large distances, g\ — *> eg 2 , 
are essentially nonperturbative effects since A ~ e -const./g e 



V. FLUCTUATING DIPOLE MOMENTS 

In previous Sections we considered the properties of the "rigid" dipoles which are char- 
acterized by a fixed absolute value of the magnetic dipole moment //. In a realistic case 
the dipole is realized as a monopole-antimonopole bound state with a fluctuating absolute 
value of the dipole moment. In this Section we generalize our approach to the case of the 
fluctuating dipole moments \i — g m r. It is convenient to describe the fluctuations of the 
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dipole moments as the fluctuations of the dipole size r at a fixed magnetic charge g m of the 
monopole constituents. In a general case the distance r can be characterized by a distribu- 
tion function F(r) normalized to unity, J +oc drF(r) = 1. We assume that the dipole gas is 



sufficiently dilute (cf. eq. (|2.1| )): < po > 3 J drr F(r) <C 1. 

The statistical sum of the gas of pointlike dipoles with the fluctuating moments is: 

oo /-N N +?° „ „ ( i N 



2 = E ^ Il( J dr i F (n) J d 3 x, J d 3 n.) exp -- £ V(g. 



N=0 ' 8=1 n -2 1 k a, 6=1 

,! nf=l 



m 1 a Tlai 9m fh ^6) ^<a> -^6/ 



(5.1) 



In the case of a fixed dipole moment, F{r) = 5{r — d), the above sum is reduced to the 
partition function (|2.2| ) with /i = g m d. 

Analogously to Section |lj we introduce the field \i perform summation over the dipole 
positions x a , and integrate out the dipole directions n. The partition function (|5.1|) is then 
represented as follows: 

1 ,* ,2 



Z = Jb X exp{- J d 3 x [-(d X ) 2 -^(v^dxl)}} 



where 



v a (f)= [ drF(r)^. (5.2) 
J r t 

o J 

Note that the action for the field \ has an additional integration over the dipole size param- 
eter r in comparison with eqs.( [2.5| , pi)| ). 

As a simple example let us consider the following distribution function, shown in Fig. 
0(a): 

F(r) = ^ 2 exp{- r -}, (5.3) 

which has a maximum at r = d, and the mean value of the dipole size, r, is f = 2d. 

Function ( |5.2j ) for distribution (|5.3|) has the following form, v a (f) = 1+d 1 2 p , and the 



expansion parameter Aa is now given by eq.( |2.7p with jj, = r g m = 2dg m : 



2„2 
in 



Afi = 16ii(d g 
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In the leading order in Aa the classical inter-particle energy can be calculated similarly to 
the previous Section: 

E^' n (R) = 4ttC J d 2 x w(x, R) = 2ir 3 (d 2 G(^) , (5.4) 

1 + d 2 r^(x 1 ,x 3 ,R) 



where T c is defined in eq. (|4.3| ) 



The leading behavior of the energy E^*' at large and small distances is given by the 
following formulae, respectively: 



7T 2 



E^ n (R) = 27r%dqR= —pdqR, #<r, (5.6) 
E^f(R) = 167T W log(-) , R » r , 



where the dipole density in the leading order in Afl is given in eq. (|2.8| ). As in the case of the 
ordinary dipoles the presence of the dipole gas leads to the appearance of the linear term in 
the inter-particle potential. The coefficient in front of the linear term is: 

o* = qof, af = 2ir 3 (d. 

At large distances the renormalization of the electric charge takes place: 

gl^etgl, e fl = 1 + ~A fl + 0{\\) . (5.7) 

The profile and contour plots for the (normalized) energy density w, eq.(|5.5|) of the dipole 
cloud are shown in Fig. § for various charge separations R. Although at small distances the 
profiles show some similarity with a stringlike object the width of such a "string" grows as 
R s t r ~ \fd~R~ for R d. Indeed, the effective width R str of the string profile at the center 
of the string can be defined as follows: 

oo 

/ x 2 3 w(0,x 3 ,R)dx 3 j 2 „ 

Rl, = = RVR , + 1M • (5.8) 

/ w(0,x 3 ,R)dx 3 



12 



Since the energy of the "string" is proportional to the volume it occupies the "string" could 
merely be responsible for R 2 tr (R) ■ R ~ R 2 (and not for ~ R) term in the potential. Thus 
the linear term in the potential has nothing to do with the observed string formation. 

It is instructive to study the behavior of the (normalized) energy of the test particle 
pair G, eq.( |5.4[ ). The plot of this function in Fig. 0(b) clearly shows that the slope of the 
function G is linear up to the distances^ R ~ 4d while the density profile, Fig. |](d), shows 
no presence of the string at R = Ad. This is another evidence in favor of the fact that the 
linear term is not caused by the string formation mechanism. 

VI. CONCLUSIONS 

We studied the static potential between oppositely charged electric particles in the point- 
like magnetic dipole gas. We found that the short distance potential contains a linear piece 
due to overlapping of the dipole clouds gathered near electric sources because of the induced 
magnetic field. This effect can not be explained by a short string formation contrary to the 
3D compact electrodynamics || and 4D (dual) Abelian Higgs model ||. The coefficient 
in front of the linear term is proportional to the density of the dipoles and to the electric 
charge of the test particles. 

As the distance between the test particles increases the potential becomes of a Coulomb 
type. At large distances the dipole effects lead to the electric charge renormalization. 

These results may have interesting applications for physics of gauge theories both in 
high and low temperature regimes. One of the physically interesting candidates of such 
theories is the electroweak model in a low temperature phase in which the formation of the 
monopole-antimonopole pairs has been observed [[L4j]. The dipole effects may also induce a 
short distance linear potential |T5| between electric charges in a high temperature phase of 
the compact U(l) lattice gauge theory in three space-time dimensions. 



3 At R ~ 4d the deviation of the classical energy from the linear behavior ( |5.6| ) is around 15%. 
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FIGURES 



FIG. 1. The Wilson loop for the static charge and anti-charge separated by the distance R. 
FIG. 2. (a) The distribution function F, eq.([y|), for the dipole parameter r, and (b) corre- 



sponding (normalized) energy G(x), x = R/d, of the charge-anticharge pair, eqs.(5.4). 



FIG. 3. The normalized energy density ( |5.E| ) of the dipole cloud around the electric sources 
in the plane (2:1,2:3). The profile and contour plots are shown for various charge separations: (a) 
R = 0.5d, (b) R = l.Od, (c) R = 2.0d, (d) R = 4.0<i. The distances are shown in units of d. 
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FIG. |(b) 
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FIG. |(b) 



21 




■3-2-10 1 2 3 



FIG. |(d) 
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